Pseudo-potential of a power-law decaying interaction in two-dimensional systems 
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We analytically derive the general pseudo-potential operator of an arbitrary isotropic interaction 
for particles confined in two-dimensional (2D) systems, using the frame work developed by Huang 
and Yang for 3D scattering. We also analytically derive the low energy dependence of the scattering 
phase-shift for an arbitrary interaction with a power-law decaying tail, V2d(p) tx p~ a (for a > 2). 
We apply our results to the 2D dipolar gases (a = 3) as an example, calculating the momentum 
and dipole moment dependence of the pseudo-potential for both s- and p-wave scattering channels 
if the two scattering particles are in the same 2D layer. Results for the .s-wave scattering between 
particles in two different (parallel) layers are also investigated. Our results can be directly applied 
to the systems of dipolar atoms and/or polar molecules in a general 2D geometry. 



Introduction: Low-dimensional strongly correlated sys- 
tems have been one of the most important subjects in 
condensed matter physics in the last few decades. From 
the many-body point of view, the standard meanfield ap- 
proximation for the 3D system is usually broken down by 
thermal fluctuation at finite temperature, while from the 
two-body point of view, the widely used first Born ap- 
proximation for a 3D weak potential totally fails in lower 
dimensional systems As a result, a proper effective 
interaction (or called pseudo-potential) becomes essential 
to go beyond the weak interaction limit of 3D scattering 
or in the long-energy limit of 2D and ID system. For the 
2D system that we want to concentrate in this paper, 
an important progress was made by Schick 2] for the 
studying of 2D bosons with a hard-disk potential, and 
the relevant applications to cold atom systems are also 
investigated by several groups [1|3@|- Recently Kanji- 
lal and Blume Q further derived a general form of the 
pseudo-potentials for all angular momentum channels of 
a short-ranged interaction. The pseudo-potential for a 
low-ranged dipolar interaction in 2D systems was also 
studied in the s-wave channel by one of us [8], but its 
extension to higher angular momentum channels is still 
unexplored, even though several important results have 
been carried out in 3D systems recently @ . 

In this paper, we systematically generalize earlier re- 
sults and apply to the systems of 2D dipolar gases: (1) 
First, we apply Huang and Yang's theory to de- 

rive a general form of the pseudo-potential in 2D sys- 
tems. Our results can be shown equivalent to Kanji- 
lal and Blume's results [IQ, derived from another ap- 
proach. (2) For a general power-law decaying interaction, 
V2d(p) P~ a (a > 2), we further analytically calculate 
the low energy dependence of the scattering phase shift, 
upto a single non-universal parameter to be determined 
by the short-ranged details of V2D- (3) Finally we ap- 
ply our results to the study of s- and p-wave scattering 
channels of dipolar interaction (a = 3), and numerically 
evaluate the non-universal parameter for a model interac- 
tion. The s-wave scattering for the two scattering parti- 
cles confined in two different (parallel) 2D layers are also 



investigated, showing a Feshbach-like resonance even at 
zero dipole moment limit. Our results can therefore be 
applied to the many-body physics of magnetic dipolar 
atoms [HI , cold polar molecules [l3| , or indirect excitons 
in a semi-conductor based double- well system [14j. 

General pseudo-potential for 2D scattering: We start 
from solving the two-particle scattering problem of the 
following 2D Schrodinger equation with total energy E: 

- ^i^(r ± ) + V 2D (p)i>(r ± ) = Etf>(r ± ), (1) 

where = j;-§^p-§^ + j?~thj? m cylindrical coordinate, 
p, is the reduced mass, and tp(rj_) is the scattered wave- 
function in the relative coordinate, = (x,y). V2d(p) 
is the effective 2D interaction, obtained by integrating 
out the transverse degree of freedom (z), and is assumed 
to be isotropic about the z-axis here. Note that we also 
have assumed that the transverse confinement potential 
is so strong that no confinement-induced resonance Q 
has to be considered here. 

Since V2d(p) is assumed to be isotropic and decays 
faster than p~ 2 in large p, the wavefunction, ip(r), can be 
always expanded by noninteracting eigenstates in large 
distance: ^ foj) = £ ~ =0 u m (k, p) £ CT=± C^Qe*™*, 

where k = y/2pE/h 2 , and u m (k,p) = A m (k)J m (kp) + 
B m (k)N m (kp) is the radial wavefunction with 
Jm(x)/N m (x) being the Bessel function of the 
first/second kind. Here A m (k), B m (k), and C^(k) 
are coefficients to be determined by boundary condi- 
tions. Similar to the 3D case [lof , we first investigate the 
short-distance behavior of such noninteracting solution 
in the leading order terms: 

u (k,p) ~ A (k) + ^-ln(^-), (2) 
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Here we define /?„ 



-~t+H m /2 with 7 _ 0.57722 being 



the Euler's constant and H m = J2k=i being the Har- 
monic number [l5[ (here Hq = 0). We note that the third 
term in the right hand side of Eq. © , resulted from the 
irregular solution, N rn (kp), is of the same order (upto a 
logarithmic function) as the first term if A m {k) ~ B rn (k). 
For a typical short-range interaction, however, this term 
can be neglected because B(k)/A(k) oc k 2m in the long 
wavelength limit. Here, since we want to derive a general 
pseudo-potential for a power-law decaying potential (see 
below) at a small (but finite) scattering energy, we will 
still keep this term for the most general application. Such 
hybridization between the regular and irregular solutions 
of noninteracting partial waves does not exist in the 3D 
case • All other terms can be shown irrelevant to the 
derivation of the pseudo-potential below. 

To derive the proper pseudo-potential, we have to ap- 
ply the noninteracting Hamiltonian on the asymptotic 
wavefunction above [lOllllI ]. and integrating over a small 
spherical area of radius p by using Green' theorem [Til ] . 
Separating contributions from the s-wave and the non-s- 
wave parts, we obtain 

~ [V 2 + k 2 } ^(r x ) = ~S(t x ) [4B (k) 



■larrup 
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(4) 



The next step is to rewrite the right hand side to be a 
function of V m (k) = B m (k)/A m (k), which is the only 
quantity related to the phase shift of the mth partial 
wave, 5 m (k). (In fact, V m (k) = — tan<5 m (fc)). In order 
to get an expression of A m (k), we have to take certain 
derivatives on the wavefunction and let p — > 0, as in the 
3D case [10|. After some simple calculation we obtain 

d r u (k,p) 
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p m u m (k,p)} ,(6) 



where we have used the following identity: 

£^(x 2m 1n(x/b)) = (2m)l(H 2m + hx(x/b)) 

As a result, by combining Eqs. (g]), (O and ([5]), we find 



^-V 2 t;(r.) >^ V,„r(r_) =Er(r_). 



E 
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where the pseudo-potential operator, V m , is 



V m = 5(t X ) 



4(m!) 2 /(2m)! 



2/i VUk)- 1 + | [H 2m + Hkp/2f3 m )] 
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Eqs. JZl)-© can be interpreted as the effect equation 
of Eq. ([T]) with the same boundary condition at origin 
(rj_ = 0) in the low energy limit (E, kp — ► 0), and hence 
is the 2D version of Huang and Yang's result in Ref. 
(To| (see Eq. (12) therein). We note that one can show 
that the pseudo-potentials derived above are equivalent 
to results of earlier work both in the s-wave channel 0, 
and in the higher angular momentum channels 




Pseudo-potential for p 01 interaction: After deriving 



the most general form of pseudo-potential for 2D scat- 
tering, we further derive an analytical closed form of the 
momentum dependence of T > m (k) for a power-law de- 
caying potential: Vzd{p) ~ U/p a as p — > oo. Here U 
measures the strength of interaction, and a > 2 is the 
decay exponent. We start from the zero energy scat- 
tering (E = k = 0) of Eq. |T]) and the radial wave- 
function, u m (0,p) — u m (p), can be calculated analyti- 
cally: u m (p) = A m Im. f||r) + B m Km. (f^ with A m 

and B m being the coefficients to be determined by the 

short-distance behavior of V2d(p)- Here £ = a/2 — 1, 

i 

A Q = (^r-) 2( , and I m (x) / K m (x) is the modified Bessel 
function of the first/second kind. In the limit of long 
distance (or weak interaction, A|/^p 5 <C 1), we have 
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which should be also reproducible by taking the zero en- 
ergy limit (k — ► 0) of Eqs. © and © (the last term of 
Eq. ([3]) can be neglected in this limit). Therefore the 
relationship between V m (k) and P rn = B m /A m can be 
easily derived to be 



V (k) 
V m {k) 



Ctt/2 



Pn 



ln(2/3oO-eMfcA a /2/3 ) 



(10) 



-27r(2f)- 



(*A a /2) 



2 m 



to! (to — 



i)! r (f + i)r(^) 



where P TO = B m j A m is the only non- universal param- 
eter, depending on the detailed shape of V2d(p) in the 
short-distance regime. Note that above results apply 
only in the low energy and/or weak interaction limit, i.e. 
k < p- 1 < £ 1/S / A a- 

For s-wave scattering channel, we can define an effec- 

-i -=i 

tive scattering length, a a = A q (2/3q£) 5 e5P ° > so that 



VoV(rx) = 



-2ttH 2 



2^1n(fca Q /2/3 ) 



<5(r_ L )V>(r J 



(12) 



if we assume the scattered wavefunction can be approxi- 
mated by a smooth function at origin after cross-graining 
the short-ranged fluctuation (for example, the meanfield 
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FIG. 1: Results for Case A: (a) and the inset are respectively 
the value of and (^o°') _1 as a function of A3. Here 
we set W = 0.1 /im (see the text), (b) and (c) show the 



calculated Vq°^ (k) and "Pj uj (k) respectively for different values 
of incident wavevector, k. 
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condensate wavefunction of bosonic particles). The re- 
sulting pseudo-potential above becomes the same as a 
hard-disk potential with an effective "radius", a a . 

The justification of the 2D pseudo-potential depends on 
the interaction strength, i.e. when ka a ~ fcA Q <C 1. 

Numerical results for the dipolar gases: In the rest of 
this paper, we will concentrate on a physical example, 
say systems of polar molecules, for the results of dipolar 
interaction (a = 3). An external electric field is assumed 
to applied perpendicular to the layer plane, inducing a 
field-dependent dipole moment, D. We consider three 
cases of scattering here: Case A: s-wave scattering be- 
tween identical bosons in the same layer, Case B: p-wave 
scattering between identical fermions in the same layer, 
and Case C: s-wave scattering between identical bosons 
or fermions in two parallel layers with layer separation d. 



The last case can be directly applied to the systems of 
multi-layer structure made by ID optical lattice In 

the rest of the paper, we will use V r 2 ( D )/(1) to denote the 
bare intra-/inter-layer interaction, with the superscript, 
to identify all the quantities obtained by either 
of them. For the convenience of numerical calculation, 
we further approximate the effective 2D interaction by 
the following analytic form: for the intra-layer interac- 

we use Vjd (p) — -pr for p > W and = ^ for 
W, where W should be about the same order of 
the layer width and is fixed to be 0.1 /im in the follow- 
ing calculation. We note that different choices of the 
cut-off, W , can bring only minor quantitative difference 
in the results of phase shift (not shown here), because 
the intra-layer interaction, V^j (p) , is assumed to be re- 
pulsive for all p, and hence no Feshbach resonance type 
resonance should be expected. For the inter-layer inter- 



tion 

P < 



where the effect of 



,. T/ {1)/ \ D 2 (p 2 -2d 2 ) 

action, we use \%'{p) = {p }+ d 2 )5 /2 , 
finite layer width is expected to be smaller since W <C d 
in a deep optical lattice. As a result, the only length 
scale associate with the our repsent model interaction is 
A3 = MD 2 /h 2 (also denoted to be ad in the literature 
|12||). For a typical molecule, say SrO, the fully polarized 
dipole moment can be D — 8.9 Debye, leading to A 3 
as large as 123.2 pm. However, for magnetic atoms like 
52 Cr, the maximum value of A3 is just about 1.03 nm. 

In Fig. [TJ we show the numerical results for Case A and 
Case B together by evaluating the original two-particle 
Schrodinger equation of Eq. ([1]): In (a) 



(0) 



we show P and 



-,(0) 


function of dipolar strength, A3. One can see 
that when A 3 is small, (Pq ^)' 1 can be quiet large, lead- 
ing to a very small scattering length, 03°' for A 3 < 0.05 

pm. However, for larger A3, 0,3 becomes proportional 
to A 3 , which is the only relevant length scale in this 
regime (i.e. the short-ranged details of the dipolar in- 
teraction becomes negligible). In Fig. [TJb), we show 

the calculated strength of pseudo-potential, V < ^ > \k), for 
different values of incident wavevectors, k. One can see 
that for a given k, (k) decreases to zero logarithmi- 
cally as A 3 — > 0, while it has a resonance- like behavior at 
a certain value of A3. Such "resonance- like" behavior oc- 
curs as So(k) = 7r/2, indicating that the interaction is so 
strong to push the wavefront of the scattered wavefunc- 
tion well-ahead of the noninteracting one. It is therefore 
nothing to do with the Feshbach resonance in 3D, and 
only results for ka^ ~ fcA3 <C 1 are correct for true low 
energy scattering. In (b) we show the results for Case B: 

as a function of dipolar strength. However, un- 
like the boson case, where the typical incident wavevector 
is determined by the condensate (i.e. system) size at low 
temperature, the typical incident wavevector for fermions 
at low temperature should be about the Fermi wavevec- 
tor (i.e. the inverse of inter-particle distance) due to the 
Pauli exclusion principle. Therefore we calculate results 
for a much larger k in (c), but similar interaction depen- 
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A 3 /d 

FIG. 2: Results for Case C: (a) is the value of a 3 /d as a 
function of As/d. In the inset: the upper one shows (P^ j ) -1 
for zero energy scattering, and the lower one shows the bound 
state energy in unit of recoil energy, Er. Note that all the 
length scales are in unit of the inter-layer separation, d. (b) 
shows the value of Vq (fc) for different incident wavevectors, 
k. Inset: the magnified plot for the first resonance. 



dence is still observed. 

In Fig. [2j we show the results for Case C: the s-wave 
scattering between particles in two different layers. In (a) 
and its inset, we show the calculated scattering length, 
03 , and the associated Pq 1 ^ as a function of A3. It 
is interesting to see that, different from the intra-layer 
case, (-Pq ) _1 diverges to negative infinity and 03 also 
diverges in the regime of small A3. Such divergence orig- 
inates from the fact that our dipolar interaction can al- 
ways sustain an inter-layer bound state in 2D system, 
even when the interaction strength is infinitely small. 
The calculated bound state energy (also in the inset) 
shows a logarithmically small binding energy for the first 
bound state, while the second bound state appears near 
A^/d ~ 71. When considering the finite size effect, i.e. 
ka is bounded below, the first resonance will occur at a 
finite dipole moment as shown in (b) (also see Ref. [|[). 
The existence of an inter-layer bound state can lead to 
a strong modification of the pseudo-potential strength 
(similar to the Feshbach resonance) , leading to some ex- 
otic many-body phases as predicted in Refs. [HHH]. 

In summary, we analytically derive the general form of 
the pseudo-potential for an arbitrary short-ranged and 
isotropic interaction in a uniform 2D system. The en- 
ergy and interaction dependence of the pseudo-potential 
is also derived analytically for an arbitrary power-law in- 
teraction. Numerical results are provided for the dipolar 
interaction, and therefore can be applied in the study of 
the 2D quantum dipolar gases. 
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